Iranian Journal of Numerical Analysis and Optimization 
Vol. 9, No. 2, (2019), pp 77-102 Poxso 


DOI:10.22067 /ijnao.v9i2.72109 yu 


Research Article 


Chebyshev pseudo-spectral method for 
optimal control problem of Burgers’ 
equation 


F. Mohammadizadeh*, H. A. Tehrani and M. H. Noori Skandari 


Abstract 


In this study, an indirect method is proposed based on the Chebyshev 
pseudo-spectral method for solving optimal control problems governed by 
Burgers’ equation. Pseudo-spectral methods are one of the most accurate 
methods for solving nonlinear continuous-time problems, specially optimal 
control problems. By using optimality conditions, the original optimal control 
problem is first reduced to a system of partial differential equations with 
boundary conditions. Control and state functions are then approximated by 
interpolating polynomials. The convergence is analyzed, and some numerical 
examples are solved to show the efficiency and capability of the method. 
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1 Introduction 


Time-dependent partial differential equations (PDEs) play an important role 
in various fields of application, such as fluid dynamics [4,28], electro magnet- 
ics [17,24], or heat transfers [10]. Often, these differential equations consist of 


*Corresponding author 
Received 17 April 2018; revised 10 February 2019; accepted 1 March 2019 

F. Mohammadizadeh 

Faculty of Mathematical Sciences, Shahrood University of Technology, Shahrood, Iran. 
e-mail: f:mohamadi093@gmail.com 


H. A. Tehrani 
Faculty of Mathematical Sciences, Shahrood University of Technology, Shahrood, Iran. 
e-mail: hahsani@shahroodut.ac.ir 


M. H. Noori Skandari 
Faculty of Mathematical Sciences, Shahrood University of Technology, Shahrood, Iran. 
e-mail: math.noori@yahoo.com 


77 


78 F. Mohammadizadeh, H. A. Tehrani and M. H. Noori Skandari 


nonlinear terms, which are challenging to solve numerically. Furthermore, in 
engineering applications, we often interest in an optimal solution of the con- 
sidered PDE with a certain objective function. This leads to the mathemati- 
cal field of PDE constrained optimization where a cost function is minimized 
and the PDE is considered as a constraint. 


Optimal control of viscous Burgers’ equation is one of the most important 
PDEs constraint optimization, which is taken into consideration and several 
papers have recently been presented in its numerical solution. Kobayashi [12] 
has concerned with adaptive stabilization and regulator design for a viscous 
Burgers’ equation by nonlinear boundary control. Yilmaz and Karasozen [29] 
by using the high level modeling and simulation package COMSOL Multi- 
physics solved the optimal control of unsteady Burgers’ equation without 
constraints and with control constraints. Also, they applied an all-at-once 
method for the optimal control of the unsteady Burgers’ equation [30]. More- 
over, in other work, they [11] transformed the optimality system for boundary 
controlled unsteady Burgers’ equation after linearization into a biharmonic 
equation in the space-time domain. Smaoui et al. [23] have dealt with the 
sliding mode control (SMC) of the forced generalized Burgers’ equation via 
the Karhunen—Loeve (K-L) Galerkin method. Hashemi and Werner [8] pre- 
sented a finite difference scheme for the one-dimensional viscous Burgers’ 
equation, which boundary conditions are taken as control inputs. Zeng and 
Zhang [31] designed a new preconditioning technique along with MINRES 
(minimal residual) in nonstandard inner product for the linear system of 
equations arising for optimal control of the unsteady Burgers’ equation. Ku- 
cuk and Sadek [13] analyzed the dynamics of the forced Burgers’ equation 
subject to Dirichlet boundary conditions by using the boundary control with 
the objective of minimizing the distance between the final state function and 
target profile along with the energy of the control. Also, they applied a 
robust technique for solving optimal control of coupled Burgers’ equation; 
see [21]. Noack and Walther [16] used adjoint techniques for efficient eval- 
uations of the gradient of the objective in gradient based on optimization 
algorithms. Marburger and Pinnau [15] showed the convergence of method 
to solve optimal control problems, where the constraints have been discretized 
by a particle method. Allahverdi et al. [2] discussed the efficiency of various 
numerical methods for the inverse design of the Burgers’ equation. An aug- 
mented Lagrangian-SQP technique depending upon second-order sufficient 
optimality condition is analyzed in [27]. In [5], a comparison of three differ- 
ent numerical methods for optimal control of Burgers’ equation is carried out. 
In [25], a Lagrangian-Newton-SQP method is presented for the solution of 
optimal control of Burgers’ equation, where the control is restricted by point- 
wise lower and upper bounds. Proper orthogonal decomposition method is 
utilized in [14] to solve optimal control problems of the Burgers’ equation. 
Furthermore, some controllability results for viscous Burgers’ equation with 
distributed controls are presented in [6]. 
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In some papers, at first the problem is discretized on one of the variables 
and then the Runge-Kutta method is applied to discrete the optimality con- 
ditions. To achieve good results, we have to take the number of points of 
discretization big enough that it can be too time consuming and has compu- 
tational complexity. 

The pseudo-spectral (PS) methods are one of the most accurate methods 
to solve continuous-time problems including ordinary differential equations 
(ODEs) and PDEs. Using the direct Chebyshev PS(CPS) method for optimal 
control problems governed by PDEs does not usually give an explicit solu- 
tion. Hence, we use an indirect CPS method for solving necessary optimal 
conditions in an optimal control problem of Burgers’ equation and we obtain 
an approximate optimal solution. It is the first time that the CPS method is 
applied by discretization of variables synchronic. By numerical examples, it 
can be seen that the CPS method is more effective than other methods and 
we can achieve the better results for the solution of optimal control problem 
of Burgers’ equation. In fact, the error of CPS method is less than that of 
other methods. 

The paper is organized as follows: In Section 2, we introduce the optimal 
control problem of Burgers’ equation. In Section 3, we give the optimality 
conditions for optimal control problem of Burgers’ equation. In Section 4, 
we utilize the CPS method to discretize optimality conditions. In Section 5, 
the convergence of the method is analyzed. Finally, the numerical examples 
and conclusion are given in Sections 6 and 7, respectively. 


2 Optimal control problem of Burgers’ equation 


The distributed optimal control problem for the Burgers’ equation can be 
stated as follows: 


1st pe a ft re 

Minimize Jly,u] = >| / (y(t, 2) — z(t, x))?daxdt + a if u(t, x)daxdt 
0 a 0 a 

(1) 


subject to 
y(t, x) 7% y(t, ©) yx(t, 2) io VYcx(t, x) = ®(u), (t, x) EQ= (0, 7) x [a, b], (2) 
y(t,a) = y(t,b) =0, texX=(0,7), (3) 


y(0, 2) => yo(x), rEQ= [a,b], (4) 


where y(.,.) is the state variable, u(.,.) is the control variable, a > 0 is the 
regularization parameter, v > 0 denotes the viscosity parameter, and ® is a 
given function. An usual selection of function ® is as follows: 
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B(u) = u win, _ 
0 uinQ—-Q, 


where 2 is the set of active controls; see [5, 25, 27]. 


At first, optimality conditions for problem (1)—(4) are given, then we indi- 
rectly develop the CPS method to achieve an approximate optimal solution. 


3 Optimality conditions 


Here, we summarize the process of achieving the optimality conditions for 
the optimal control problem (1)—(4). 

Let Lj,.(Q) be the set of all functions that are Lebesgue integrable on every 
compact subset of 2 and let W*:?(Q) be the linear space of all functions 
w € L?(Q) having weak derivatives Dw in L?(Q) for all multi-indices a of 
length | a |< k. Assume that C§°(Q) is the space of differentiable functions 
of every arbitrary order on 2. The closure of C§°(Q) in W*?(Q) is denoted 
by W"?(Q). Moreover, we define Hk (Q) = Wh? (Q). 

Definition 1. Let y € [j,,.(Q) and let some multi-index a be given. If 


loc 


a function w € Lj,,.(Q) satisfies 


/ y(x)D°v(x)dx = (—1)!«! i w(a)u(x)dx, for all vu € CH°(Q), 
2 2 


then w is called the weak derivative of y (associated with a). 


Let H = L?(Q) and V = Hj(Q) be Hilbert spaces. We make use of the 
following Hilbert space 


W(0,T) ={¢ € L7((0,T];V) : de € L7((0, 7]; V*)}, 


where V* denotes the dual space of V. The inner product in the Hilbert 
space V is given with the natural inner product in H as 


(¢,v)v = (Ov )x, for OW € V, 


where ¢’ and w’ are the derivative of ¢ and w, respectively. 


Definition 2. Every function y € W(0,T) that satisfies 


< y(t), d >vev +v(u(), ov + (¥Oye(t), da = (fF + ®(u))(t), Da 
gEV, te [0,T], (y(0),x)x =(yo,x), xe AF, 


where 
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<ult)orviv= [ node f rf.o'de, ne V. 
Q Q 
is called a weak solution for Burgers’ equation (2) and conditions (3) and (4). 


In order to achieve the optimality conditions, the operator e : X —> Y 
is introduced by 


e(y, u) = (er(y, u), e2(y, u)) = (yt — Vex + YYe — f - &(u), y(0) — Yo): 


where y(0) = y(0,2) , X = W(V) x L?(Q), and Y = L?(V) x H is identified 
with Y* = L?(V*) x H the dual of Y. We use L?(V) for L?(X;V). Now, 
the optimal control problem (1)—(4) can be transformed as the following 
minimization problem with equality constraints 


minimize J(y,u) 


subject to e(y,u) = 0. 


Theorem 1. Let (y*,u*) be an optimal solution of (1)-(4). Then there exist 
Lagrange multipliers p* : [0,T] x Q —> R and r* satisfying the first-order 
necessary optimality conditions 

E' (y* ,u*,p*, A*) =0, e(y*,u*) =0 
with the Lagrangian 


Ly, U, DP, A) = J(y,u) — (ex(y, u), P) £2(V*),L2(V) = (e2(y, u), A)A- 


Hence, first-order optimality conditions lead to the following optimality 
system (see [9, 26, 27]) 


Yt —VYce + Yo = Pu), (t,z) € Q, 

Mts + yPe=ya-y, (4a)€ Q, 

y(t, a) = y(t, b) = 0, te &, 

y(0, x) = yo, xe Q, (5) 
p(t, a) = p(t, b) = 0, te &, 

p(T, x) = 0, re Q, 

au+p=0, (t,x) € Q. 


From the last equation of system (5), we have 


u= —*p. (6) 


By using (6), we can express (5) as follows: 
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Yt —VYze + YYe = O(-Sp), (2) € Q, 


Pit VPrx + YP = Ya — Y; (,2)€ Q, 

y(t,a) = y(t,b) =0, te &, i 
p(t, a) = p(t, b) = 0, te X, 

y(0, 2) = yo(x), re Q, 

p(T, x) = 0, re QQ. 


Now, we utilize the CPS method to solve optimality conditions (7) and obtain 
an approximate optimal solution for the optimal control problem (1)—(4). 


4 CPS method to approximate the optimal solution 


Chebyshev polynomials are in center of the CPS method, and hence we first 
introduce these polynomials and then demonstrate the CPS method to solve 
optimality conditions (7). 

Chebyshev polynomials [1,18] are orthogonal polynomials, which play an 
important role in the theory of approximation. These polynomials are defined 
on [—1, 1] as follows: 


These polynomials can be generated from the following recurrence relation- 
ship: 


Tj+1(2) = 2xT; (x) = T;-1(2), j= 1,z2e€ [-1, 1). (9) 
It is possible to give an explicit expression of Chebyshev polynomials as 
T (a) = cos(7 cosa). (10) 


It is now easy to verify that Ty(.) has N zeros within the interval (—1, 1) 
as Zp = Cos ses k = 1,2,...,N. However, in this paper, we use the 
CGL points on [—1,1] which are as follows: 

TI : 
=— — 0<j<N 11 
Lj cos 5, j ‘ (11) 
where they are the roots of (1 — x”)T\(x). For interpolating in the CPS 
method, the following Lagrange polynomials are utilized: 


N 


N 
G— Ly 2 1 
Ly (x2) = J — T; (vz) Tj (x), k=0,1,...,N, x € [-1,1], 
ON xe [-1,] 
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where 


_ s4, 7 =0,N, 
aa 1, 1<j<N-1, 
and we have 
j=k, 
0, j#k. 


To use the CPS method, the variables of system (7) are transformed to in- 
terval [—1, 1] by the following linear transformations: 


(12) 


t==t+ ot t€[0,T], t€ [-1,1], 
ba ipa ; (13) 
C= E+ , «© €[a,b], & € [-1,1). 
2 2 
By (13), the optimality conditions (7) can be written as follows: 
Y; = W1 YEeyee Ds Vale, &), Yer(E, i) 3 (t, z) € [—1, 1] x [-1,1 ’ 
P= pe (OE @), PUG) IG, @)s Peat, Z))-s (t, Z) € [—1, 1] x [=1,1 ) 
Y(#,-1) = Y(#,1) =0, te [-1,1], 
P@,—1) = P(é,1) =0, te [=l, 1), 
Y(-1,£) = Yo(2), ze {-1,1], 
P(1,z) =0, te [-1,1], 
(14) 
where 


P(t,Z) = t+ z 
(4,2) =P(St+ 51> ae 
= T. T b-a b+a 
Yi(t, Z) = t T 
dine) yal 5 te D x 2 13 
b-a b+a 
y, tT) = 
0(Z) yo( 9 9 Ms 
T 2 2 1 
Y, Yes) = = (| (——)?vYae z 
tO PY You) =F (GaP — ays — 8-4), 
T 2 49 
be (Y, P, Pz, Pra) = ( )°uPre YP;+Ya—-—Y 
2 b-a b— 


We assume that w,; and v2 have bounded and continuous derivatives with 
respect to their arguments. Hence, there exist constants M, and M2 such 
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that 


lea (Y, P, Ye, Yes) — Ya (% P, 2, Yen) |< My (Iv - ¥|4 (P= Pl), 
(15) 


Ivo (¥, P, Ps, Paz) — 2 (YP, Ps, Pes) |< Mo (IY -¥1+|P- PI). 


Now, to approximate the optimal solution, we utilize the following polyno- 
mials interpolating: 


AG, #) = Die 0 Lie —9 a ayy L,(t)L;( ms (16) 
PN(,z) = oes 9 Oly Li (E)L; (2), 


where ayy and by ; 


(12), we have 


for 7,7 = 0,1,...,N, are the unknown coefficients. By 


Yj =ah 
5) aa (17) 
P (ti, Z;) = bi 
where : ; 
a _ Tj 
i, = —cos(—), = —cos(—). 
cos(—), 2; cos( =") 


To express the derivatives Y;¥(.,.), Y(.,.), YAX(,.), PX (..) PH(.,.), and 


PX(.,.), we can use the matrix multiplication D = (D;;) and get 
N(F = NaN 
Ye" (tp, Ek) = dic Uk aij,Dpi, 
YA i 0) = j= (0) ap; Drj: 
N(Fom) — 
Yaz (tp, Tx) >= dia 10) ap Day, (18) 
PN (tp, Ex) =~ yy 9 bin Dpis 
Px (tp, Ex) = Ljco boy Das 
Pah (tp, Be) = Dj 0 be Dig: 
where 
Lk . i 
el = Fh 
Mj | tr t 
tk 0<j=k<N-1 
_ = _ #2? = = | 
Dey = Tile) =) Sua ih, (19) 
a JF k = 0, 
2N2+4+1 
— j=k=N, 


and D = D- D = (Dx) where Dy = eg DesDy, ky § = 0,1,...,N. In 
fact, multiplying by the matrix D transforms a vector of the state variables 
at the CGL points to the vector of approximate derivatives at these points. 
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Now, by relations (16), (17), and (18), conditions (7) can be written as the 
following discrete form: 


N 
pone 0 aj), D pi — th (a Qk» ae pee ) GY Digs 4 0 a’ Dui) =0, 
= ND tbe (ad ‘pk? epee 0b Phir N= o BN Drs) = 0, 
Gy = 0, boy = boy =0, GGL, = Yo(Zx), ON, =0,k,p =0,1,...,.N. 
(20) 


an 


By solving the system of algebraic equations (20), we can obtain pointwise 
and continuous approximate optimal solutions as (16) and (17), respectively. 
Also, by (6), the approximate optimal control can be given as 


UN (6,2) = = LW OL@, Ga) € 1,1) 41), 


Moreover, the optimal value of functional (1) can be approximated as follows: 


N ON 
T = /o b-a_ b+a _ 
= 5 S wattle — = (6p +1), 5 Bet ) acy (21) 


k=0 p=0 


where a a 1oN., and ws,s =0,1,...,N, are the quadrature weights of the 


numerical approximation (21). For even N, the weights are 


1 
Wo = WN = Fe 
a a 2 = (22) 
2 —- 
ty = WN-2 = H Fo Taz 008(“*); 8 = 1 QyivagS, 
and for odd N, 
wo = Wn = ws 
—_ cal 2 
i a 1 2798 N-1 (23) 
— — 2 = as 
Ws =WN-s= WH Xj20 TC Tage c0s( Sy SN 2 ang 


The double prime in the weights formula denotes the first and the last ele- 
ments have to be halved. 


Remark 1. Notice that the approximation solutions (16) can be written 
as the following Kronecker (tensor) product form: 


N WN 
Y(é,2)~Y% (2) =) ay Li (OL;(Z) = LE, 2)Ty 
i=0 ss 


P(é,@) » P’ (é,2) Sy (OL; (#) = L(E,z)Tp, 


1=0 7=0 
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L(t, £) = [Lo(t)Lo(@)--- Lo) Ln (@)L1 (6) Lo(&) --- Li) Ln (Z)--- 
Ln (é)Lo()--- L(t) Ln (2)| 
= [Lo(t) Li(é)--- Lv (#)] ® [Lo(#) L1(Z) --- Ln(2)] 
= L(t) ® L(z), 


and the symbol “@” denotes the Kronecker product. Also 


-N ..=N -N..=N ~N ~N 
Ty = [€99°+- 19° Ay + GN GNI’, 
7N ..7N ZN ...3N 7N 7N 
Tp = [boo ++ bon big ++ Oty + b+ byl 
where the symbol “’ ” denotes the transpose. All the components of Ty and 


Tp are unknown and our aim is to obtain them (A similar representation can 
be seen in [22, 32, 33]). 


5 The convergence and error analysis 


In this section, first we give the following definition and then analyze the 
convergence of the presented method. 

Assume that © = [-1,1] x [-1, 1]. We apply C"(Q) for the space of func- 
tions with continuous derivatives of rth order. 


Definition 3. [19] The function W : R* > Rt with the following prop- 
erties is called a modulus of continuity if 


i) W is increasing, 


) 
ii) lim, 49 W(z) = 0, 

ili) W(z1 + 22) < W(21) + W(2z9), for any 2 and z2 € R*, 

iv) there exists a constant c such that cW(z) > z for allO << z <2. 
Some important modulus of continuity can be defined as 


W(z)=2%, O0<a<l. (24) 


Now, assume that B? is the unit circle in R?. We say that a continuous 
function f(-,-) on Q admits W(.) as a modulus of continuity, if the following 
value is finite 
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Suppose that space Cjj,(B) includes all functions f(-,-) on B? with contin- 
uous first order partial derivatives, and it equips with the following norm: 


FCs Dlaw = MFC Dlhoo + WFEGs Dlheo + IfeGs )lleo + [feCs lw + Ife, lw: 


(26) 
Also, we define the space Cjj,(Q) as follows: 
Co a = {f(.,)¢€ C'(Q) : for all (t,@) € ©, there exists a map¢ : B? > Q, 
. (#,%) € Int(o(B*)) and f o o(-,-) € Cy(B’)}. (27) 


It can be proved that if Q = (i. , Int(¢;(B?)) for some ¢1,...,¢7, then 


f(,:) € Ch Ww (Q) if and only if f o 4;(-,-) € Cw(B?) for each i = 1,...,1. 
Moreover, ch Ww(@) with norm 


lfC daw = Sisoat )\li,w (28) 


is a Banach space (for more details, see [19]). Now, we show the space of all 
polynomials of total degree at most 2N on 2 by Pol(N, N,Q), that is, 


Pol(N, N,Q) = {n(t, Z) = Yow Z : (£,%) €O, yi; € R}. 


Theorem 2. For any f(-,-) € Cjy(Q), there is a polynomial n(-,-) € 
Pol(N, N,Q) such that 


Coci 


If) — mG Mlle $ Fe 


—Ww( (29) 


ay 


where c, = ||f(-,-)lla,w and co is a constant that depends on W(-), but inde- 
pendent of N. 


Proof. The proof is a result of Theorem 2.1 in [19]. 


Now, to guarantee the existence of solution for the system of algebraic 
equations (20), we relax it as follows: 


N 
Lie 0 ajy.D pi — tr (aN nk ce Sel 0 cae ee Sel =0 aN Daj, je 0 aN Dus) | 


“VE ~ 
< YY wigs), vk =1,2,...,.N-1 


ae 
Dixo bi Dpi — pe (a.,b it wh Dri dss 0 DG, o BY Das) | 


< Wass aha) p,k=1,2,...,N—1 
oe = AAW ats gna): |arwl < = ys STEEL ae 
b50 pol < OW ads gNai)> [on wl < SAW (sh syri)> pP=9,1,...,N, 


ag, — Yo(Zx)| s Me Wop 3N=1)> lbvnl < = MEW onl 3N=1)> k=0,1,...,N, 


(30) 
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where N is sufficiently big and W(-) is a given modulus of continuity. Since 
a W (ath) = 0, any solution (,5)~ = ((a%,,8%,);p,k = 0,1,...,N) 
for system of algebraic inequalities (30) is a solution for system of algebraic 
equations (20) when N tends to infinity. 

Now, we show that system (30) is feasible, that is, it has at least one 
solution (G, b) yn. 
Theorem 3. Let (Y(-,-), P(-,-)) be a solution for system (14), where Y(-,-) 


and P(-,-) are in Cy,(Q). Then there is a positive integer K such that, for 
any N > K, the system (30) has a solution as 


(@, b)n = (Gon, Opn )s Dp, & =0,1,...,.N). (31) 


Moreover, (a,b) Nn satisfies 


Y (tp, Be) — Gp < sea W (sy=z)) p,k =0,1,...,N, 
|P (Gp, Bx) — BA < sia W (sx): p,k=0,1,...,N, 


where L is a positive constant independent of N. 


Proof. Assume that 7,(-,-) and m2(-,-) in Pol(N — 1, N,Q) are the best 
polynomial approximations of Y;(-,-) and P;(-,-), respectively. By Theorem 
2, we get 


1 ie ia et asta a (7,2) 


(33) 
IP(E2) — mE B)lleo < g@2zW og), (3) € 0 


where 7; and 72 are two constants independent of N. We define 


V(é,@) =Y(-1,2) + fi, m(7,2)dr, (62) € an 
P(E,2) = P(-1,2) + f°, m(7,2)dr, (2) <Q, 
and 7 7 7 
GY th), by SH Pete. PHO 1M (35) 


We show that (a,b) = ((@),6),); pk = 0,1,...,.N) satisfies system (30). 
By (33), (34), and (35), for (£,%) € Q, we get 


T 


|Y(é2) —Y(é,2)| = icc ©) —m(r,#))dr| < [. [Ye(7, 2) — m(r,®)| dr 


V1 1 & 211 1 
<= < : 
< aay) | ars an 1 (gn —7? 


(36) 


Also, by a similar procedure, for (¢,%) € Q, we gain 
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PE,8) — BG a)| < 55? Wg). (37) 


Now, by relation (34), the functions Y(-,%) and P(.,Z), for any % € [—1, ll, 
are polynomials of total degree at most 2N. Hence, their derivatives at CGL 
nodes fo, f,...,£y are exactly equal to the value of polynomial at the nodes 
multiplied by the differential matrix D, defined by (19). Thus, we have 


N 
So ain Dpi = lps te) So biLDpi = tn Bas P; k= 0, 1, acars /N. (38) 


Therefore, by relations (15) and (36), we get 


N N N 
N —N -N Ff 

5 Ai,Dpi — V1 ADK» Oe 2 a; Dj; s a; Dk; 

i=0 =0 j=0 


< ale 
+ |¥eltp; Tx) ais Soa Pas Dom Dy 
= 71 (tp, Lx) = Yeltasite)| 
a V1 (Y (bos a) Pls Ba) Vo (Fp, he), Yor bo Be)) 
a N N 
<1 | apis pis Day Pras d | nj Drs 
j=0 j=0 
< |n1 (tp, Zk) — Ye(tp, Ze)| + Mi (|Y (tp, Fe) — GY] + |P (tp, Be) — BPI) 
V1 i 271 1 272 1 
< Ww M 
~ 2N-1 (soi) t (ea) ta SG) 
= 1 (2M, +1) + 2Mi7 1 _ 
—_ IN —1 Woy) p,k=1,...,N 1, (39) 


where M, and M2 are Lipschitz constants of 1, and W2, respectively. More- 
over, by a similar process, we obtain 


N N N 
Di Dpi — Ya | Ape, Boks D | Bey Daar D | bps Drs (40) 
i=0 j=0 j=0 
(2M + 1) + 2Mo27V1 1 = 
< aN 1 Woy) p,k=1,...,N—-1. 
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Further, for boundary conditions, we get 


IA 


7 ie- Yo(x) P31, g1~V (Hie) eV eleo= Yay) 


271 1 


< Ws), B= 0,1. (41) 
Also, for all p,& = 0,1,...,N, we have 

P(1,ax)| = |PU.ax) — PUL.a8)| < se WG) 42) 
ad a a Oi Oak 

#@,-1)| = |? &,-1) -¥G@,-1)| < 271 wr) _) (43) 
- ie ec = ~2N-1 ‘2N-1” 

VE | = V(&,1) -VE vy] < 271 _ yr} _) AA) 
ee P= ON 1 SON —1" 

~ ~ 272 1 

P ,-1)| =|PG,-1 Pf, -1)] < (4 
,-0)| = |PG 1) - P| < Wigs), C48) 

BS, us Boe oo = 272 1 

PG, 1)| = |P, 1) - PGs < Ww 4G 
| =|PG.1) - PG.D| < 55 W gy) ) 


Hence, if we select K such that 


max{71(2M, + 1) + 2My72,72(2M2 + 1) + 2Mo71, 21, 272} < VN, 


for N > K, then by (39)-(46), pair (@, b)y satisfies system (30). 


Now, we show that the sequence of solutions of problem (30) and sequence 
of their interpolating polynomials converge to the solution of the problem 
(14). 


Theorem 4. Let {(a,b)y = (CARA p,k = 0,1,...,N}S Lx be a se- 
quence of solution of system (30) and let {(YN(-,-), PX (-,-))}%_¢ be their 
interpolating polynomials sequence defined by (16). Also, we assume that, for 
any © in [—1,—1], the sequence {(Y" (-1, #), P% (-1,2), eae -), PRG, Wer 
has a subsequence {(YNi(—1, x), PN‘(—-1, 2), eG, -), Pee, -))} S25 that uni- 
formly converges to (f° (€), 63° (Z), A1¢, +), Ae, :)) where A1(-,-), Ae(-.+) € 
C72(O), %(-), O3(-) € C2([-1,1]) and limjo0 Nj 

= oo. Then the pair 


(¥(¢,2), P@) = (Jim Y" (Ea), lim P™,2)), (47) 
for (t,Z) € Q, is a solution of system (14). 
Proof. By attention to the assumptions, we have 


¥ (2) = OP @) + fo, u(r B)dr, (48) 
P tz Te 


(E,@) = $3°(@) + fi, A2(7, B)dr. 
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We first show that (Y(£, 2), P(E, Z)) for f € [— 
satisfy system (14). Assume that (Y(-, Zz), 
does not satisfy the first equation of (14). 


that 


1, Yand@ = az, k =0,1,...,N 
P(-,%,)) for some k = 1,..., N—1 
Then, there is a 7 in (—1,1) such 


¥e(r, Bx) — da (¥ (7, 2e), Plr, Bx), For, Bn), Voor, Bx) # 0. 


Since CGL nodes {tp}}9 when N -+ oo are dense in [—1,1], there is a 
sequence {ity, 2, such that 0 < ly, < Nj and limo tin, =rT. Thus 


On the other hand, since lim;_,,, VNi Wisc—a 


IN; 1 ) = 0, by (30), we obtain 


a 


in (Yilaiy, 2 k) — i (¥ (tiy, »k); Pity, , Bx), Ya(tiy, » Fe), Yen (fix, sFx))) = 0, 
which is a contradiction to (49). Thus (Y(£,Z), P(é,#)) (for all  € [—1, 1] 
and = Z,%, k = 1,...,N —1) satisfies the first equation of (14). By a 
similar procedure, we can show that it satisfies the second equation. Also, 
it can be easily proved that (Y(-,%x),P(-,£x)), for k = 0,1,...,.N, satisfies 
the boundary conditions. For example we show that Y(—1, %,) = Yo(%,) for 
k=0,1,...,N. We have 


0 < |¥(—1, Z,) — Yo(Zx)| =| lim Y**(—1, Z%) — Yo(Zx)| 
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= lim |Y%*(—-1, %) — Yo(Zx)| = lim |a}it — Yo(Ze)| 
1 CO 1—0o 


VNi 1 
ws 2N,—-1 


IA 


lim 
izoo 2N;—-1 


y=0. 


Hence Y (1, Zz) = Yo(%,) for all k = 0,1,...,N. Now, we know that 
nodes {Zk }iLo when N — ov, are dense in [—1,1]. Therefore the pair 
(Y(-,-), P(-,:)) defined by (47) is a solution for (14) on Q = [-1,1] x 
[—1, 1]. 


6 Numerical examples 


In the following examples, we use the Levenberg—Marquardt method (a quasi- 
Newton method) for FSOLVE command in MATLAB software to solve alge- 
braic system (20). 
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Example 1. Consider problems (1)-(4), where T = 1, a = 1, v = 0.01 


and 0.05, yo = sin(472), ®(u) = u, and z(t,x) = 0. The approximate opti- 
mal value of objective function computed by the CPS and LPS [20] methods 
for vy = 0.01 and 0.05 and N = 10, 20, 30 and 40 are shown in Table 1. 
We observe that our numerical results are better than the results of the LPS 
method. In Figures 1 and 2, we show the obtained approximate optimal state 
and control for v = 0.05 and different values of N. 


(a) The approximate optimal state for N = 10 and v = 0.05 (b) The approximate optimal state for N = 20 and v = 0.05 


YNitx) 


(c) The approximate optimal state for N = 30 and v = 0.05 (d) The approximate optimal state for N = 40 and v = 0.05 


Figure 1: The approximate optimal state for Example 1 
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(a) The approximate optimal control for N = 10 and v = 0.05 
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(b) The approximate optimal control for N = 20 and v = 0.05 


(c) The approximate optimal control for N = 30 and v = 0.05 


(d) The approximate optimal control for N = 40 and v = 0.05 


Figure 2: The approximate optimal control for Example 1 


Table 1: Comparison of objective function values for Example 1 


vy =0.01 


vy =0.01 


v = 0.05 


vy = 0.05 


Presented method 


LPS method [20] 


Presented method 


LPS method [20] 


0.033014881174862 


0.0828638100277 


0.016911085761598 


0.01590006952876 


0.040440356851832 


0.0620867108909 


0.013730446517693 


0.01519309308076 


0.029728232559725 


0.0466282421253 


0.015082505388501 


0.01519227846689 


0.029005091596013 


0.0463124455511 


0.015073940981453 


0.01519176630695 


Example 2. Consider optimal control problem (1)—(4). Let T = 1, a = 0.05, 
vy = 0.01, and z(t,2) = yo(a) where 


yo(x) = 


1 


, «é€ (0,5), 


0 otherwise. 
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is the initial condition and 


u, won (0,T) x (4,4); 


0, otherwise. 


Table 2 shows the approximate optimal values of objective function for dif- 
ferent values of N by the CPS method. In Figures 3 and 4, the optimal state 
and optimal control are presented, respectively. 


(a) The approximate optimal state for N= 15  (b) The approximate optimal state for N = 20 


(d) The approximate optimal state for N = 40 


(c) The approximate optimal state for N = 30 


Figure 3: The approximate optimal state for Example 2 
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(a) The approximate optimal control for N = 15 (b) The approximate optimal control for N = 20 


(c) The approximate optimal control for N = 30 (d) The approximate optimal control for N = 40 


Figure 4: The approximate optimal control for Example 2 


Table 2: Approximate values of objective functional for Example 2 


N Jn 

20 | 0.106416588965259 
25 | 0.085739365971895 
30 | 0.083129631634715 


Example 3. Consider optimal control problem (1)—(4). Assume that T = 
10, a= 1, ®(u) = u, and z(t,x2) = yo(a) where yo is defined in Example 2. 
We apply the CPS method to this problem. In Table 3, the values of cost 
functional using CPS method for different values of N are listed. In Figures 
5 and 6, the approximate optimal state and optimal control are illustrated, 
respectively. 
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Table 3: Approximate values of objective functional for Example 3 


N Jn 

10 | 0.795502231705757 
20 | 0.732649373291908 
30 | 0.605003194783743 


(a) The approximate optimal state for N = 10 (b) The approximate optimal state for N = 20 


(c) The approximate optimal state for N = 30 (d) The approximate optimal state for N = 40 


Figure 5: The approximate optimal state for Example 3 
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(a) The approximate optimal control for N =10  (b) The approximate optimal control for N = 20 


UNtx) 


(c) The approximate optimal control for N = 30 (d) The approximate optimal control for N = 40 


Figure 6: The approximate optimal control for Example 3 


Example 4. Consider the optimal control of Burgers’ equation (1)—(4) with 
T = 10, a = 0.1, v = 0.01, ®(u) = u, and the desired state z(t,x) = 
Yyo(x) where yo = exp(—2) sin(27x), a € [0,1]. The numerical results are 
displayed in Table 4 for different values of N by using our method. In Figures 
7 and 8, we show the approximate optimal state and control, respectively. 


Table 4: Approximate values of objective functional for Example 4 


N Jn 

20 | 0.121081851037199 
30 | 0.116980721067958 
40 | 0.115852005217238 


yNax) 


yNatx) 
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(a) The approximate optimal state for N = 10 (b) The approximate optimal state for N = 20 


yN(t.x) 


(c) The approximate optimal state for N = 30  (d) The approximate optimal state for N = 40 


Figure 7: The approximate optimal state for Example 4 
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UN(tx) 


(a) The approximate optimal control for N= 10 _(b) The approximate optimal control for N = 20 


UNtx) 
UNit.x) 


(c) The approximate optimal control for N = 30 (d) The approximate optimal control for N = 40 


Figure 8: The approximate optimal control for Example 4 


7 Conclusions 


In this paper, we proposed an efficient Chebyshev pseudo spectral method 
to solve the optimal control problem governed by Burgers’ equation. By 
applying this method, we discretized the optimality conditions and obtained 
a system of algebraic equations. We achieved a good approximate optimal 
solution with good accuracy. We analyzed the feasibility and convergence of 
the presented method. 
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